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Abstract 

A systematic method for the computation of finite temperature (T) crossover 
functions near quantum critical points close to, or above, their upper-critical 
dimension is devised. We describe the physics of the various regions in the 
T and critical tuning parameter (t) plane. The quantum critical point is at 
T = 0, t = 0, and in many cases there is a line of finite temperature transi- 
tions at T = T c (t), t < with T c (0) = 0. For the relativistic, n-component 
(j) 4 continuum quantum field theory (which describes lattice quantum rotor 
{n > 2) and transverse field Ising (n = 1) models) the upper critical dimen- 
sion is d = 3, and for d < 3, e = 3 — d is the control parameter over the entire 
phase diagram. In the region \T — T c (t)\ <C T c (t), we obtain an e expansion 
for coupling constants which then are input as arguments of known classical, 
tricritical, crossover functions. In the high T region of the continuum the- 
ory, an expansion in integer powers of y/e, modulo powers of lne, holds for 
all thermodynamic observables, static correlators, and dynamic properties at 
all Matsubara frequencies; for the imaginary part of correlators at real fre- 
quencies (lj), the perturbative y/e expansion describes quantum relaxation at 
Tuo ~ ksT or larger, but fails for Tko ~ y/eksT or smaller. An important prin- 
ciple, underlying the whole calculation, is the analyticity of all observables as 
functions of t at t = 0, for T > 0; indeed, analytic continuation in t is used to 
obtain results in a portion of the phase diagram. Our method also applies to 
a large class of other quantum critical points and their associated continuum 
quantum field theories. 
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I. INTRODUCTION 



The study of finite temperature crossovers in the vicinity of quantum phase transitions 
is a subject with a long history p^pl|, but many aspects of it remain poorly understood. 
The structure of the crossovers is especially rich for the case where the quantum critical 
point extends into a line of finite temperature phase transitions, and there is a reasonable 
qualitative understanding of all the regimes. While there have been quantitative calculations 
of crossover functions in special cases [|[||[1^,[18|-|2(J there is no complete, general theory of 



these crossovers, especially for the case when the quantum critical point is below its upper 
critical dimension. 

In this paper, we shall provide a new systematic and controlled approach to the quantita- 
tive computation of these crossover functions. Our method is quite general: it should apply 
to essentially all quantum critical points in the vicinity of, or above, their upper-critical 
dimension. 

Recently, O'Connor and Stephens [[EJ have also studied crossovers near relativistic 
quantum-critical points below their upper-critical dimension. They found it necessary to 
introduce a non-standard extension of the field-theoretic renormalization group. We will 
comment on their results (and of others) in Section fVl. 

In this paper, we will show that it is possible to devise a simple strategy, completely with 
the framework of standard field-theoretic methods, which provides a systematic computation 
of the required crossovers. We shall describe how our method can be extended to arbitrary 
orders in an expansion in powers of the interactions, but we shall only provide here explicit 
computations at low orders. One of the main virtues of our method is that it clearly separates 
contributions of fluctuations of different physical origins: critical singularities of the T = 
quantum-critical point, and those of the finite T classical phase transition, are accounted 
for at distinct stages of the calculation. 

We shall present most of our discussion in the context of a continuum quantum field 
theory (CQFT) of a n-component bosonic field <p a (a = 1 . . . n; we will drop the index a 
except where needed) with 0(n) symmetry and with the bare, imaginary time (r) action 



S = J Q drj d d xl- 



(<w + (v0r + K c +t o 



- \ d.i; 



We have set % = Hb = 1, measured length scales (x) in units in which the velocity of 
excitations c = 1, and introduced the bare mass m 2 oc + to and the bare coupling ii (the mass 
term has been separated so that the T = quantum critical point is at t = 0). This field 
theory describes the low-energy physics in the vicinity of the quantum phase transition in 
the (i-dimensional transverse-field Ising model (for n — 1) or the 0(n) quantum rotor model 
(for n > 1). The generalization of our method to other quantum field theories, like the dilute 
Bose gas, or the models for onset of antiferromagnetism in Fermi liquids, is straightforward 
and will also be discussed. 

We begin our discussion by reviewing the expecting scaling structure of S for the case 
where the quantum critical point is below its upper critical dimension. At T = 0, S describes 
the usual 4 theory in d + 1 dimensions, and its upper critical dimension is d = 3; for d < 3, 



there is an essentially complete understanding fl23"| , p4| of the critical properties of this theory 
in an expansion in powers of 
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e = 3 — d. 



(1.2) 



The definition of the renormalized theory requires a field scale renormalization Z, a coupling 
constant renormalization Z±, and a renormalization of 2 insertions in the critical theory 
Z 2 . In terms of these, we define us usual 

t = t Z/Z 2 ; (1.3) 

t is a measure of the deviation of the system from its T = quantum critical point. Precisely 
the same renormalizations are also sufficient to define a finite theory at non-zero T, even in 
the vicinity of the finite T phase transition line, as we shall explicitly see in this paper. 

We show a sketch of the phase diagram of S as a function of t and T in Fig [T] [|21[| . 
We have assumed in this figure, and throughout the remainder of the paper that the model 
S is in its scaling limit i.e. the coupling uq is at its fixed point value, and all ultraviolet 
cutoffs have been sent to infinity after an appropriate renormalization procedure. There 
is a finite temperature phase transition line at T = T c (t), and all other boundaries are 
smooth crossovers. All of the physics is contained within universal quantum-critical crossover 
functions, which we now briefly describe. We will consider the behavior of the dynamic two- 
point susceptibility, x{k,^) (k is a spatial momentum, and u is a frequency) obtained after 
analytic continuation of the susceptibility 

X {k,iuj n ) = \\otx fJ T dTe^ kx -^(4>(x,T)4>(0,0)) (1.4) 

which is evaluated at Matsubara frequency u n = 2mrT. We will consider the scaling behavior 
of x f° r t > and t < separately, and then discuss the relationship between the two cases. 
(i) t > 0: The susceptibility x obeys the scaling form || 

„ , , / he V ( k B T\ TI T fhek hu A+\ 

where we have momentarily re-inserted all factors of h, ks, and c, r] is the usual field 
anomalous dimension of the T = 0, d + 1 dimensional theory, and is a fully universal, 
complex-valued, universal scaling function. Notice that there are no arbitrary scale factors, 
and x is fully determined by two parameters, A + and A, which are properties of the T = 
theory. The first of these, A + , is the true energy gap above the ground state, while the 
second, A, is the residue of the lowest quasi-particle excitation; they obey 

A+~f ; A~tf v t>0 (1.6) 

where v is the usual correlation length exponent of the T = theory (all Greek letter 
exponents in this paper refer to those of T = the quantum-critical point, and not to 
those the finite T phase transition). We provide a computation of the values of the T = 
parameters, A + and A + , in Appendix |B[ 

The factors in front of in (|1.5|) have been chosen so that is finite at A + /k B T = 0. 
All scaling functions defined in this paper will share this property. 

We also emphasize that although the scaling ansatz (|1.5|) contains dynamic information, 
its form and content are quite different from the dynamic scaling hypotheses applied near 
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classical phase transitions fl25j. In these classical systems, a single diverging correlation 



length, £, is used to set the scale for k and u; the analog of (1.5) is then a scaling function 
of two arguments, k£ and co>£ Zc , where z c is the classical dynamic critical exponent. In 
contrast, the quantum crossover result Q1.5Q is a function of three arguments, the extra 
argument arising because the quantum critical point has two relevant perturbations (T and 
t). Further, the identification of universal scale factors, and indeed the conceptualization 
of the physics, is much more transparent when T is used as the primary energy setting the 
scale for other perturbations. Only in the immediate vicnity of the finite T phase transition 
at T c (t), \T — T c (t)\ <C T c (t), does ( |1.5| ) collapse into a scaling function of two arguments, as 
has been discussed in Refs P,|26||. 

(ii) t < 0: Now the T = ground state breaks a symmetry with 



t) = N S, 



al 



T 



0,t < 0; 



;i-7) 



here N$ ~ (—t) 13 is the condensate, which we have arbitrarily chosen to point in the a = 1 
direction, and (3 = (d — 1 + r])vj2 is the magnetization exponent of the T = theory. Now 
Nq can serve as the parameter which determined the field scale (replacing A), and we need 
an energy scale which determines the deviation of the ground state from the t = 0, T = 
quantum critical point. For n = 1, there is a gap, A_, above the ground state, which satisfies 
our requirements; we have therefore || 



N*(hc) 



A li r 1+ri (k B T) 



2- V 



hck Kuj A_ 
knT 1 knT 1 kuT 



t <Q,n 



1. 



;i.s) 



For n > 2, there is no gap above the ground state, and we use instead the spin-stiffness, 
Ps ~ (—tp d ~ 1 ^ as a measure of the deviation from the quantum critical point; in this case 
we have the scaling form 



N 2 

Xa(k,Uj) = — 
Ps 



he 



\k B T) 



d-X\v/{d-l) 



ihc) d ~ 2 Ps 




t < 0,71 > 2. 
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A computation of the T = parameters N , A_ and p s is provided in Appendix |B|. The a 
dependence in Q1.9Q accounts for the difference between fluctuations transverse and longitu- 
dinal to the condensate orientation. Again, is finite at A_/k B T = 0, or at p s /T d - 1 = 0, 
and all subsequent scaling functions will share this property. The finite temperature phase 
transition at T c (t) is contained entirely within the scaling function this transition ap- 
pears as a point of non-analyticity of as a function of A_/T or p s /T d ~ l . An immediate 
consequence is that the value of T c can be determined precisely in terms of the T = energy 
scale; we found, in an expansion in powers of e, that 
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k B \fe 
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and 
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1/2 



2n 2 (n + 2) 



0(e,e 



(l+e)/(l-e) Al{2vY 



n>2 (1.11) 



The various higher order contributions are all universal, but arise from very different physical 
effects; we will discuss their origin later in the paper. In the upper critical dimension (e = 0), 
these formulae are modified by replacing e by 1/ ln's : so for n = 1, T c ~ A_ ln 1 ^ 2 (l/A_) etc. 
For the cases n — 1, d — 1 and n > 3, d = 2 it is known that in fact T c = i.e. long-range 
order is present only at T = 0, and disappears at any non-zero T. For these cases, it is clear 
that the above results for T c , and other results obtained in the e expansion, cannot be used 
in the region labeled III in Fig [I]. However, the results of this paper can still be usefully 
applied to the remainder of the phase diagram of Fig |l|. 

Although we have used two separate scaling forms to describe the behavior for t < and 
t > 0, there is a crucial connection between them. Notice that there is no thermodynamic 
singularity at t — provided T > 0. This implies that the observable x(k,uj) (and indeed 
all other observables) must be analytic as a function of t at t = as long as T > 0. This 
principle will serve as an extremely important constraint on the calculation in this paper; 
indeed, our method is designed to ensure that analyticity holds at each order. Further, our 
results in the region t < 0, T > T c (t), were obtained by a process of analytic continuation 
from the t > 0, T > region (see Fig [!]). The ground state for t < has a spontaneously 
broken symmetry, and hence cannot be used to access the symmetric t < 0, T > T c (t) region 
in perturbation theory; instead it is more naturally accessed from the disordered side with 
t > 0. A similar procedure of analytic continuation in coupling constant was used recently 
in exact determinations of quantum critical scaling functions in d = 1 . 

Before turning to a description of our method in Section [[ A| , we highlight one of our 
new results. A particularly interesting property of CQFT's at finite T is the expected 
thermal relaxational behavior of their correlators in real time. This behavior cannot be 
characterized simply to a field-theorist who merely considers correlators of S, defined as a 
CQFT in imaginary time with a spacetime geometry R d x Si (the tensor product of infinite 
(^-dimensional flat space with a circle of circumference 1/T) . In real frequency, the thermal 
relaxational behavior is characterized by the fact that lim^o im x(^? w)/u is expected to be 
finite, with the limiting value proportional to a relaxation constant. In Ref f2~H| , the quantity 



R 



«x(o,o) 



;i.i2) 



uj=0 



was introduced as a convenient characterization of the relaxation rate. By (|L 
obey a scaling form given by 



k B T 

h 



k R T 



Tr must 



(1.13) 



for t > 0, and similarly for t < 0. In particular, in the high T limit of the CQFT PT] (this is 
region II of Fig|l|), T R is iksT/h) times \I/r+(0), which is expected to be a finite, universal 
number. Unfortunately, we shall find that the perturbative expansion discussed in this paper 
cannot be used to obtain a systematic expansion for T^. A self-consistent approach, with 
damping of intermediate states, appears necessary and will not be discussed here. In the 
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high T limit, we shall show that the non-self-consistent approach fails for frequencies of 
order y/eT/h or smaller. To avoid this difficulty, let us define an alternative characterization 
of the damping at frequencies of order ksT/h by 

(1.14) 

uj=k B T/h 

This will obey a scaling form identical to that of Tr. We found in the high T limit that 

+ 0(e 3/2 ), (1.15) 

where Li2 is the dilogarithm function defined in (|2.36| ). The result will be compared with 
an exact result for Tr? for d — 1, n — 1 in Section |IIB| . Also in the exact result at d — 1, 
n = 1 we find Trt/Tj- = 0.98170018 ... in the high T limit, and we expect a similar ratio 
close to unity in all cases in region II of Fig [[J 

The following subsection contains a description of our approach. Readers not interested 
in the details of its application to S, can read Section [1 A| , and skip ahead to Section [TV| 
where we give our unified perspective of this and earlier work. 

We will revert to setting H = ks = c = 1 in the remainder of the paper. 



t rt = -X(0, 0) Im 

UJ 



RT 



(n + 8W 



- + 7T 2 + 6Li 2 (e 
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A. The Method 



The origin of the approach we shall take here can be traced to early work by Luscher |27 
on the quantum 0(n) non-linear sigma model in 1 + 1 dimensions. Subsequently, a related 
idea was employed by Brezin and Zinn- Justin |28| and by Rudnick, Guo and Jasnow p9f in 



their study of finite-size scaling crossover functions in systems which are finite in all, or all 
but one, dimensions (also referred to as the d — * and d — > 1 crossovers). The quantum- 
critical crossovers are clearly related, but now involve d + 1 —> d. We shall show here that 
the latter problem can be successfully analyzed by essentially the same method as that used 
for the former. There are some new subtleties that arise in a limited region of the phase 
diagram, and we will discuss below how they can be dealt with. 

We will describe the method here for the special case of the action S (Eqn (|1 . 1|) ) with 
d below its upper critical dimension i.e. d < 3. The central idea is that at finite T, it is 
safe to integrate out all modes of <fr(k,u) n ) with a non-zero Matsubara frequency, u n ^ 0, to 
derive an effective action for zero frequency modes <fr(k, uj n = 0). All modes being integrated 
out are regulated in the infrared by the term in their propagator, and so the process is 
necessarily free of infrared divergences; further, the renormalizations of the T = theory 
also control the ultraviolet divergences at finite T. To be specific, let us define 

T r 1 / T 

$(x) = —= \ dr<f>(x,T), (1.16) 



and its transform in momentum space $(/c) = J d d xe thx $>(x). Then from S, we can deduce 
an effective action for $(/c) after completely integrating out the <fi(k,uj n ^ 0) (we have set 
the coupling uq at the fixed-point of its /3-function — see Appendix H): 



6 



c _ 1 



1 r rl d h ~ 
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2 J (2tt) 

+ 1_J d%dj2d% d ^ ^ ^ ^ k3)Hki)Hk2)HhM _ ki _ ^ - * s ) + . . . . (L17) 

The couplings C 2 , C4, . . ., are computed in a power series in e, with the coupling constants 
renormalized as in the T — 0, d+1 dimensional critical theory. This procedure will remove 
all the ultraviolet divergences of the quantum critical point. However, the ultraviolet di- 
vergences of the finite T, (^-dimensional 4 theory remain; fortunately these are very simple 
as the 4> theory is super-renormalizable for d < 4 |30| . In particular, d < 3, ultraviolet 



divergences are associated only with the one-loop 'tadpole' graphs; so let us define 

cm ^ m+ (^y ^ (k ,- k M,- kl) + (118) 

where the ellipses refer to "tadpole" contributions from higher order vertices like Cq, C§ .... 
Similarly, there will also be tadpole renormalizations of C4 to C4 by higher order vertices, 
and so on. These new vertices, C 2p , p integer, are now free of all ultraviolet divergences (for 
3 < d < 4 there is a second classical renormalization at the two-loop level which must be 
accounted for; we will ignore this complication here and deal with it later in the paper). 
They are also automatically free of infrared divergences as we are only integrating out modes 
with a finite frequency. Indeed, these vertices must obey the scaling forms: 

C 2p (h) = T^-^) (*±J < ^ , (1.19) 

for t > 0, with the a set of quantum scaling functions; the subscript Q emphasizes that 
these scaling functions are properties of the T = quantum-critical point, and distinguishes 
them from classical crossover functions we shall use later. Similar scaling results hold for 
t < and we will refrain from explicitly displaying them. For our subsequent discussion it 
is useful to define a set of couplings K, R, and U, which can be obtained from the C 2p , and 
which play an important role in our analysis: 

R = C 2 (0) 

dC 2 (k) 



K= 1 + 



dk 



2 



k=0/ 

U = C 4 (ki = 0) (1.20) 

It is clear that R, K, and U obey scaling forms that can be easily deduced from ( |1.19| ). 

We also note that the couplings R and C 2v {ki) are guaranteed to be analytic as a function 
of t at t — 0; this is because only finite frequency modes have been integrated out, and their 
propagator l/(k 2 + uj\ + 1) is not singular at t = in the infrared {k = 0) as u n 7^ 0. This 
analyticity will be of great use to us later. 

Assume, for the rest of this section, that we know the ^q± functions (we will provide 
explicit computations of some of them later in the paper). We are now faced with the 
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seemingly difficult problem of computing observables in the $ theory with the action «S e //- 
This is a theory in dimension d close to 3 (not d close to 4), and one might naively assume, 
that this problem is intractable. We shall now argue that, in fact, it is not. The argument 
is contained in the following two simple, but important, observations. 

(1) Consider a perturbation theory of «S e // in which the propagator is l/(Kk 2 + R), and we 
expand in the non-linearities C2 P (ki), p > 2 in powers of e. This expansion differs from an 
ordinary expansion in powers of e only in that the mass R and coupling K itself contain 
corrections in powers of e; alternatively one could also treat K — 1 as an interaction, and 
work with the propagator l/(k 2 + R), but it is essential to keep the mass R in the propagator. 
Such a procedure is guaranteed to be finite in the ultraviolet. This follows immediately from 
the statement that the renormalization of the T = theory (which we carried out while 
obtaining <S e // from S) is sufficient to remove ultraviolet divergences even at non-zero T. In 
other words, the momentum dependencies in the C^/cj) must be such that all ultraviolet 
divergences cancel out. 

(2) The action <S e // is weakly coupled over the bulk of the phase diagram in the t,T plane, 
and so the procedure in (1) leads to accurate results for physical observables. Only in the 
region \T—T c (t)\ <C T c (t) (drawn shaded in Fig|l|) is a more sophisticated analysis necessary, 
which will be described momentarily. To verify this claim, consider the values of the low 
order couplings in S e ff at t — 0, but T finite; we will find later that 

R ~ eT 2 -" 
K ~ T- r > 

U „ e T e - 2 ^ (1.21) 

for t = 0. (For the present purpose, we can neglect all C 2p for p > 2 as they are all of order 
e 2 .) A dimensionless measure of the strength of non-linearities in <S e // is 

11 ! ^M/^^fl-^l «1 fort = 0. (1.22) 



K d/2 R (4-d)/2 

The above dimensionless ratio is simply that appearing in the familiar Ginzburg crite- 



rion 



|31~f . So a simple perturbative calculation is adequate for t = 0. For t > 0, the behavior 
of perturbation theory can only improve as the mass R becomes larger, which decreases 
the value of the above dimensionless ratio; as a result the perturbative calculation describes 
the crossover between the quantum-critical and quantum-disordered regimes of Fig 0. For 
t < the perturbation theory is initially adequate, but eventually becomes unreliable in the 
region \T - T c (t)\ < T e (t). 

As the above results contain some the key points which allowed the computations of this 
paper, it is useful to reiterate them. We describe the nature of our expansion, excluding the 
shaded region \T — T c (t)\ T c (t) of Fig The first step is to obtain an expansion for the 
'mass' R of the uo n = mode: we outlined above a procedure which yields a series in integer 
powers of e. Then, generate an expansion for the physical observable of interest, temporarily 
treating R as a fixed constant independent of e; this will again be a series in integer powers 
of e, but strong infrared fluctuations of the u n = mode lead to a singular dependence of 
the latter series on R. Finally, insert the former series for R into the latter series for the 



8 



physical observable. In the quantum-disordered region (Fig [T]), R is of order unity, and the 
final result remains a series in integer powers of e. However, in the quantum-critical region, 
R is of order e (Eqn (|1.21|) ) and the result, by ( |1.22j ), is a series in integer powers of y/l, 
with possibly a finite number of powers of lne multiplying the terms. It is important to 
note that the final series in powers of y/e is obtained from the original series in powers of e 
only by a local re-arrangement of terms i. e. given all the terms up to a certain order in the 
e series, we can obtain all terms below a related order in the \fe series. 

Finally, let us turn our attention to the troublesome region |T — T c (t)\ <C T c (t). We 
expect this region to be dominated by the classical fluctuations characteristic of the finite 
temperature transition, and hence to be well described by the following action Sc, which is 
a truncated form of S e ff\ 



1 



d d x 



^(v$) 2 +-$ 2 +^y 

2 v 2 4! 



1.23) 



The couplings above were defined in ( |1.20| ). We have implicitly performed tadpole renor- 
malizations where necessary to remove ultraviolet divergences of the classical theory. An 
immediate consequence of the super-renormalizability of the classical Sc is that all observ- 
ables are universal functions of the "bare" coupling constants K, R and U. So for example, 
we have for the static susceptibility 

1 1 / Kh 2 TTJ \ 

X(k^ n = 0) = -(Mk)\ 2 ) = j^c Kd/2R w /2 ) (1-24) 

where is a universal crossover function with no arbitrary scale factors. In fact, the 
crossover function ^c has been considered earlier in Ref where it was dubbed the 
tricritical crossover function for entirely different physical reasons (we emphasize that this 
terminology is purely accidental - we are not dealing with any tricritical point here). The 
computation of tricritical crossover functions is a logically separate problem from those 
considered in this paper, and we shall have relatively little to say about them here. We 
shall simply treat them as known, previously computed functions; for completeness, we 
tabulate some results on these functions in Appendix |A]. Notice that the arguments of the 
classical crossover function ^c i n ( |1-24| ) are themselves quantum-critical crossover functions, 



as follows from (|1.19p and ( |1.20|) . Indeed, inserting ( |1.19|) and ( 1.20|) into ( |1.24|) , we get a 



scaling form completely consistent with (|1 . 5| ) . The critical temperature, T c , is determined 
by the condition x(0, 0) = oo; in general, this is not equivalent to the requirement R = 
(although this does turn out to be the case at the one- loop level), but is instead given 
by the point where the scaling function ^c diverges. This condition leads to equation 
TU I \K d l 2 R( i ~ d ^ 2 ) = constant where the constant is determined by the point where ^c 
diverges as a function of its second argument; it is the solution of this equation which leads 
to the (^(e^+^/t 1 ^)) corrections to the result for T c reported in (|1.10|) . 



To summarize, in the region \T — T c {t)\ <^ T c (t), the physics is described by universal 
crossover functions, which are "crossover functions (^c) of crossover functions (\1/q)". 
The ^>q functions are properties of the quantum-critical point, and it is the burden of 
this paper to compute them; these functions then serve as arguments of known classical, 
tricritical crossover functions (^c)- 
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Finally, we note that for d just below 3, it is also necessary to include a coupling V§ 6 
in S to get the correct infrared behavior; we have ignored this complication for simplicity; 
moreover, as V ~ e 3 and this effect is present only at a rather high order. 

The outline of the remainder of the paper is as follows. In Section [II] we will compute 
properties of S for d < 3. The discussion is divided into a one-loop computation of static ob- 
servables in Section [II A| and a two loop computation of dynamic observables in Section [II B 



Section |TJ will present a general discussion of the properties of models above their upper 
critical dimension: the modifications necessary in the scaling forms and the explicit compu- 
tation of crossover functions. Finally Section [TV] will review the main results, discuss their 
relationship to previous works, and point out directions for future work. A number of details 
of the calculations appear in the appendices. The tricritical crossover functions appearing 
in ( |1.24| ) are in Appendix |A|. In Appendix JTJ we compute the T = parameters that appear 



in the scaling forms. Details of the finite T two-loop computations of various quantities are 
in Appendix |C[ 

II. CROSSOVER FUNCTIONS OF S BELOW THREE DIMENSIONS 

A number of crossover functions for the model S were introduced in Section |. We give 
formal expressions, valid to two-loop order, for all of these quantities in Appendix y. In 
this section we will evaluate these expressions and show that they obey the required scaling 
forms order-by-order in e. We will discuss the behavior for general values of t only to one- 
loop order in Section [II A| . We will limit our explicit two-loop results to a few important 



quantities at the critical coupling t — 0; these will appear in Section [UB| . 

The same basic trick will be repeatedly used to evaluate the frequency summations in 
Appendix 0: we will always subtract from the summation of a function of u n , the integration 
over frequency of precisely the same function. The resulting difference will always turn out 
to be strongly convergent as a function of momentum in all d. So, for example, we have: 



TJ2 



d d q 1 f de f d d q 



(2n) d e 2 n + q 2 + a 2 J 2tt J {27f) d e 2 + q 2 + a 2 

d d q 1 f d d+1 p 



TJ2 



(27r) d e 2 n + q 2 + a 2 J (2vr) d p 2 + a 2 
d d q 1 1 
(2n) d e yft*+^/T~[ 



(2.1) 



Notice that the integrand on the right hand side of the last equation falls off exponentially 
for large q, and the integral is therefore convergent for all d. 



A. One loop results 

We will begin (Section [II A 1|) by determining the coupling constants R, K and U of S c i 
which are the arguments of the tricritical crossover functions. Subsequently, we will present 
results for observables of S: the susceptibility x (Section |II A 2|) and the response of the 
system to a field that couples to the conserved 0{n) charge (Section [11 A3| ). 
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1. Coupling constants of S c i 



We begin by using ( |C7| ) and ( |1.20|) to obtain an expression for R, valid to one-loop order: 



R = t Q + u 



n + 2 
6 



d d q 



T 



(27TY \ e % el + q 2 + t g 2 



d d+i p i 
(27r) d +V 



(2.2) 



We have set Z = 1 at this order, and will implicitly do so in the remainder of this section. 
We now apply the identity (|2.1|) , perform the momentum integrals over p, express in terms 
of the renormalized t using ( |1.3| ). Finally we express Uq is terms of a renormalized coupling 
g defined by 



Z 2 

g = /i S d+1 — m , 
Z4 



(2.3) 



where \x is a renormalization momentum scale, = 2/(T(d/2)(Air) d ^ 2 ) is a phase space 
factor, and the values of the renormalization constants are tabulated in (|B1|). This gives for 
t > 0: 



, n + 2 
/> i 1 - -^-<) : /''// 



/ n+2 , t 
t[l + e— In — 

V 2 n + 8 /i 2 



n + 2 
12 

+ 2 



r ( 2 — - ) 1 



n + 2 
6 



T 



-1 + -) t 1 ^/ 2 
t 



J^2 



+ eT 2 



n + 2 
n + 8 



J^2 



In the second equation, we have evaluated at the fixed point value g = g* (Eqn 
then expanded to order e. The function Fa is given by 



F d (y) 



d d k 



1 



1 



1 



S d+1 J (2n) d [V^Ty e V^-l k^ + y k 2 
Let us also note here the values 

(4vr) 1 / 2 r((ci+l)/2)r(d- l)C(d- 1) 2tt 2 



(2.4) 
, and 

(2.5) 



F d (0) 



T(d/2) 



F 3 (0) 



(2.6) 



It is now easy to see, using (|B6|) that the result ( |2.4|) for R can be written in the scaling 
form 



R 



T 2 ~vAl {A 2 



2^2 



A 



+ e 



n + 2 
n + 8 



J>2 



(2.7) 



This result is consistent with the scaling postulated in ( |1.19 ). At this order, the exponent 
7] = 0, and verifying the powers of f] in front requires a higher order computation. 

We now wish to extend this result for R to t < by analytic continuation from the t > 
result. First, we need to verify that the t > result for R is analytic at t = 0. To do this, 
we first rewrite (12. 41) in the form 
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R = t\l 



n + 2, T 
e m — 

n + 8 /i 



+ 2 

+ eT 2 

n + 8 



t 

2^2 



where 



G + (y) 



ylny 



+ 4 



k dk 



y/k 2 + y eV^+w 



+ 27T^ 



(2.8) 



(2.9) 



The first term in (]2.8|) is clearly a smooth function of t; if we can now show that G + (y) 
is a smooth function of y near ?/ = 0, we will have established the analyticity of R at 
t = 0. Performing an integration by parts of the integral in (|2.9| ), followed by an elementary 
re-arrangement of terms, we can manipulate the result for G + (y) into the following form: 



G + (y) 



ykvy 
2 

ylny 



+ 27Ty/y-A J dkln 

+ 2>KyJy + 2 J™ dk (^Jk 2 + y - k - In 
sinh( v / Fl r y/2) 



k 2 + y 



y 



2VW+1 



-4 



cifc In 



2 4VF+1 



(2.10) 



The first integral can be done analytically, and we find that all the potentially singular terms 
cancel. Our final expression for G + (y), valid for y > is: 



G+(y) = |-4 / o °°rffc (in 



k 
2 



4^ 



(2.11) 



It should now be evident that ( |2.11 ) is a smooth function of y at y = 0; the integrand 



involves only even powers of ^k 2 + y, and its integral is a smooth function of y. Indeed, it 
is not difficult to explicitly extend the above result to z — —y < 0. Divide the integral into 
the regions k < y/z and k > y/z; the integrand remains unchanged in the second region, 
while in the first region the sinh function becomes a sin function — this gives us the function 
G-(z) as the analytic continuation of G + (y) to z = —y < 0: 



v ; 2 



z sinh 1 y/z — 4 



dk In 



o 



k 



sin(y/z - k 2 /2) 



4 J™dk fin 



k 



sfz~^¥/2 
smh(y/k 2 - z/2) 

y/W^/2 



k z 
~2 + Ay/W 



(2.12) 



We can now combine the above results, to obtain an expression for R valid for both signs 
of t, and which is smooth at t — 0: 



R = t[l + e 



n + 2T\ 

In - + eT 

n + 8 n) 



,n + 2 



G — 

n + 8 \T 2 



t 



G( y ) = 9(y)G + (y)+9(-y)G4-y) 



(2.13) 



A plot of the function G(y) is shown in Fig 0. As expected, the plot is smooth at y = 0. 
However, the alert reader will notice that there is in fact a logarithmic singularity in G (R) 
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at y = —2tt (t = — 47r 2 T 2 ) where the argument of the logarithm in (|2.12 ) can first change 
sign. However, this singularity is of no physical consequence as it occurs when the system 
is already in the ordered phase (Fig [I]), and the above expressions can no longer be used; 
the transition to the ordered phase happens when t ~ — eT 2 . More precisely, we see from 
the tricritical function in Appendix [A] that the value of T c (t) is determined by the condition 
R = 0; applying this to (gig ) and using (|BT3|) , (|BTT| ), (pT7|) we get the results ([LlOD , 
( |1 . 1 1| ) for T c . We will discuss the physical significance of the limiting behavior of R and G 
in various regimes in Section |II A 2 . 



Next, we turn to the computation of U. First, we obtain from ( |C5| ) and ( |1.2(J| ) the 
expression 



U = u Q - u 



n 



6 



^E 



d d q 



1 



(27TY (e 2 + g 2 + to) 2 



(2.14) 



As in the computation of R, we can set Z — 1, r\ — in the one-loop approximation. Now 
note that 



T E 



d d q 



(2ny (e 2 + g 2 + t y 



d d q 



1 



(2-nY (e 2 + g 2 + t ) 2 
d d q 1 



T 



(2ny (g 2 + 1 ) 2 
d d g 1 



( 27r )<m ( p 2 + t ) 2 , 
1 



dy \J (2ir) d y/q^Ty e y/^+y _ 



1, 



where ?/ = t^/T 2 . Using the definition of the function G + in ( |2.11|) , we get finally 
d d q 1 



(2.15) 



T E 



(2tt)<* (e 2 + g 2 + t Q f 



T~ e S d 



+i 



J^2 



0(e] 



[2.16) 



Note that we have analytically continued G' + to G' and obtained an expression which is 
manifestly analytic at t = 0. Inserting ( p. 16 ) into ( 2.14j ), expressing wo and t in terms of the 
renormalized g (Eqn ( |2.3| )) and t (Eqn ( |1.3[ )); this yields 



U = ifg [1 + 
Evaluating at g = g* (Eqn 

u 6<r 



n + 
~67 



+ 



n + 8\ \i 2e g 2 



6e 



-- + G" 

e 



y2 



(2.17) 



and expanding to order e 2 , we get finally 



n + 8 



1 + e 



'3(3n + 14) I s 
. (^ + 8) 2 "2. 



+ eG" 



t 

J^2 



(2.18) 



Note that the /z dependence has dropped out at this order, and this result is consistent with 
the scaling form ( |1.19|) and ( |1.20| ). 

Finally, it is clear that the coupling K — 1 at one loop. 
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2. Susceptibility 



The one-loop susceptibility follows immediately from the result (|C9"|): 
X -\k, ico n ) = k 2 + cu 2 n + R-e ZkTVR 



(2.19) 



As R is analytic at t = 0, so is this result for \- It is also clear that this result obeys the 
scaling forms ( |1.5| ), (|1.8p , and ( |1.9| ), given that the result (|2.7| ) for R obeys ( |1.19| ). 

The result (|2.19| ) gives us a prediction for the T and t dependence of the correlation 
length £: 



r 2 = /?- e ('^W/v7? 

n + 8 



(2.20) 



By examining the limiting behavior of £, we can obtain a physical interpretation of the 
regimes of the CQFT associated with the t = 0, T = quantum-critical point, as shown in 
Fig 

(I) T <C A + : low T limit of CQFT; paramagnetic phase 



From (J272y) and (gj) we find 



r 



Ai + e 



n + 2 



n 



T(87rTA + ) 1 / 2 e- A+ / T + . . . 



(2.21) 



So the correlation length, and the physics, is dominated by its t > 0, T = behavior, with 
exponentially small corrections due to a dilute number of thermally excited quasiparticles. 
(II) T > \t\ zv : high T limit of CQFT 

In this case, the leading behavior of £ from fl2.20| ) and ( |2.7f ) is 



ri + 2\ 2tt 2 T 2 
n + 8 



'6e(n + 2) 
n + 8 



1/2' 



(2.22) 



The scale of £, and indeed of all the physics, is now set by T. The ratio £ _2 /T 2 is a universal 
number, obtained above for small e. The reasons for the appearance of the a/c terms were 
discussed earlier in Section [1 A\ as also noted there, notice that there were no such terms in 
the region I. This series for £ is not useful as it stands, as it has the unphysical feature of 
changing sign for physically interesting values of e and n. 

As noted earlier, to this order in e, the phase boundary T = T c (t) in Fig [l] is determined 
by the condition R = and yields the values for T c given in fll.lOp and ( Jl.llp (the order 
e l l 2v corrections follow from assuming the scaling form (|2.7|) and the fact that F% has an 
expansion in integer powers of e). The result for % and R in this subsection are not valid in 
the region |T — T c (t) \ T c (t), where, instead, we have to insert the results for R, K, and U 
in Section [II A 1| into the tricritical crossover functions of Appendix |A]. In the ordered phase 
we have a second low T limit of the CQFT (region III of Fig |1|) where again the properties 
are dominated by scales set by the t < 0, T = ground state (p s for n > 1, and A_ for 
n = 1). A separate analysis with a spontaneously broken symmetry is necessary here; it can 
be easily performed by our methods, but we have not presented it in this paper. 
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Let us now turn to dynamic properties. At this order, the self energy has no momentum 
or frequency dependence; as a result the imaginary part of the susceptibility contains only 
delta functions at real frequencies: 

lm X (k, u) = ^y[% - e{k)) - S(u + e(k))] (2.23) 

with e 2 (k) = k 2 +t;~ 2 . This is clearly an artifact of the one-loop result, as the spectral density 
is required on general grounds to be non-zero at all frequencies at any non-zero temperature. 
The two-loop computation of the imaginary part of the susceptibility in Section [II B| will 
not suffer from this defect. 



3. Response to a field coupling to the conserved 0{n) charge 

The 0(n) symmetric action S possesses a set of n(n — l)/2 conserved Noether charges. 
In this subsection we will examine the susceptibility, xh, associated with an external field H 
which couples to one these charges. This analysis is motivated primarily by recent work || 
on the d = 2, 0(3) sigma model of two-dimensional quantum antiferromagnets, where 
this susceptibility is the response to an ordinary uniform magnetic field. Here, we will 
complement the earlier 1/n expansion results || by the e expansion. It is also worth noting 
here that what we have denoted here as the ordinary susceptibility x(k, to) is the staggered 
susceptibility of the quantum antiferro magnet. 

Let us orient the field H such that it causes a precession of <fi in the 1 — 2 plane. The 
time derivative term in ( |1.1| ) is then modified to 



(drfa - iHfaf + (d T <J) 2 + iH<j> x ) 2 + J2(d T <f) a ) 



a=3 



(2.24) 



The susceptibility, xh, is the second derivative of the free energy with respect to variations 
in H. We can evaluate this using the method described in Section [TS] and Appendix to 
first order in uq we obtain 



XH = 2T y f gg ? 2 + *o-4 
XH e ^J (2rrY(e 2 n + q 2 + t ) 2 



+ 2T 



d d q 



-2u 




(27r) d q 2 + R 
1 1 



d 1 ^ 4 + Qi + to ql + R 



d d q 2 qi + t c 



sn 2 n 



(2*y (n* + q 2 + t f 



(27i) d \q 2 + R 



T 



d d q 2 



d d+1 p 1 
(2vr) a! +V 



(27r) d (q 2 2 + Rf 



^2.25) 



Evaluating the frequency summations and the momentum integrals, expressing in terms of 
the dimensionless coupling g (Eqn ( |2.3|) ), and expanding some of the terms to the needed 
order in e, we obtain from ( ^.25| ) 
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Xh 

rpd-l 



t \ 2T(1 - d/2) 
Vd \T*) + (4vr) d / 2 



TV 



g fn + 2 
2 



6 



+ 9 



n 



6 



t) ® d (t 2 



" ,T 2 / T 2 e 2T 2 



(2.26) 



where the function F^y) was defined in (|2.5| ) and 

d d fc 1 



Qd(y) 



2r(i 



(2vr) d 2sinh i ( v ^ 2 T^/2) (47r) d / 2 



(2.27) 



A number of important results now follow from (|2.26 ) and ( |2.27| ); as the analysis is quite 



similar to that in Section II A 1 , we will omit the details: 
[i) After expressing in terms of the renormalized t by to — t(l + gin + 2)/(6e)) (Eqns ( |1.3| ) 



and ( |B1| )), we find that the poles in e cancel to order g. 
(ii) The resulting expression for xh is then analytic as a function of t at t — 0. This follows 
from the previously established analyticity of R at t = 0, and the fact that Qd{y) is analytic 
at y — 0. The result ( p. 26 ) can therefore be used both for t > and t < 0. 
(m) For t > 0, express £ in terms of the true energy gap A + (using (|B6p), and evaluate 
fl2.26j ) at the fixed point coupling g = g* (Eqn ( P2Q ). All dependence on the renormalization 
scale n disappears, and \h then satisfies the scaling form 



Xh 



T 



[2.28) 



where is a universal function, easily obtainable from Q2.26 ). A similar result holds for 
t < 0, where the renormalized energy scale is now the spin stiffness p s , related to t by ( |B11|) 
and (1ET71 ). 

We will be a little more explicit at the critical coupling t = 0. First, we have 



Qd{0) 
Q'd(o) 



4(d-l)r(d-l)C(d-l) 
(4vr) d / 2 r(d/2) 

+ 0(e) 



An- 



(2.29) 



Using these results, 
Xh 1 



D and Q, we get from (|2~26|) : 



- Ve 



n + 2 



1/2 



/2(n + 2) 

+ e 777 ^ + 0.310311256 

\3(n + 8) 



+ G(e 3 / 2 ) 



(2.30) 



at t = 0. At the physical value for two-dimensional antiferromagnets, n = 3, e = 1, the 
successive terms in ( |2.3U| ) oscillate in sign, and do not become smaller - so a direct evaluation 
does not yield a useful numerical estimate. 

B. Two loop results 



All computations in this subsection will be limited to the critical coupling t — 0. 
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fc=0 



Two loop results for the values of the static quantities R, x _1 (0> 0) an d dx~ l {k, 0)/dk 
are presented in Appendix 0. Our main purpose in obtaining the results is that they provide 
an explicit demonstration of the consistency of the method proposed in this paper: all 
ultraviolet and infrared divergences cancel as required, and the results take the form of a 
systematic series in powers of y/e, along with a finite number of factors of log(e). 

In this subsection we will limit our discussion to dynamic observables, in particular those 
related to lmx(k, u). Two loop contributions make an important qualitative difference in 
that the delta functions peaks in (|2.23| ) are broadened due to dissipative thermal effects. 

We begin with the expression (|C5|), retaining only the terms dependent upon the external 
frequency, setting the coupling u to its fixed point value, and keeping terms up to formal 
order e 2 : 



_ x . 2 2 2(n + 2) T 2 r 



2(n + 2)T 2 ^ f d 3 qi d 3 q 2 1 

(27r)3(2vr)3(gf + ?(e n ))(g 2 2 + 2F(fi„)) 

x — — — — — — — - + terms independent of cu n , (2.31) 

((« - ?i - <?2) 2 + o(oj n - e n - ll n )) 

where at the critical coupling t — (compare ( jClOQ ) 

a(e n ) = e 2 n + R6 tnfi (2.32) 

Notice that the finite frequency propagators in (|2.31|) have only their bare mass which 
vanishes at the critical coupling t = 0, while the zero frequency propagator has a fluctuation- 
induced mass of order eT 2 . Clearly, this distinction is an artifact of our method which treats 
the zero frequency modes in a manner distinct from the finite frequency modes. However, the 
distinction is unimportant in an expansion for physical quantities as a series in y/e (modulo 
logarithms of e), for the finite frequency modes have a minimum value of u 2 = 4nT 2 , which 
overwhelms any mass term of order eT 2 we might consider adding to their propagator. In 
other words, ( |2.31| ) provides the leading frequency- dependent contribution to x _1 (u;„) in an 
expansion in y/e for all physically allowed values of uj n . 

We are interested here in the value of lmx(k,uj) for real frequencies uj. In principle, 
this can be obtained by analytic continuation from the values of the susceptibility at the 
Matsubara frequencies. However, and this is a key point, there is no guarantee that the 
analytically continued result will also be a systematic series in y/e, valid for all values of 
uj. In fact, it is not difficult to see that the analytically continued result is valid only 
for uj ^> y/eT. This condition can be traced to the ambiguity in the mass term for the 
finite frequency propagators discussed in the previous paragraph; while this ambiguity is 
unimportant at the Matsubara frequencies, a simple estimate shows that it strongly modifies 
lmx(k,uj) for uj ~ y/eT. As a result, we are only able to obtain here systematic results for 
Im%(fc, a;) for uj y/eT. The e dependence of the very important low frequency limit uj — > 
for finite T remains an open problem. Similar difficulties were also encountered earlier in 
the 1/N expansion of the same problem, where the expansion broke down for uj ~ T. 
Here, we are able to explore the region y/eT <C w < T, and in particular have systematic 
results for uj ~ T. 

In the remainder of the section we will therefore restrict our attention to w > y/eT . 
Under these conditions we can drop the mass R even from the zero frequency propagators 
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while computing the imaginary part (all infrared divergences controlled by a finite R occur 
only in the real part). Evaluating the frequency summation in ( p. 31 ), analytically continuing 
to real frequencies, and taking the imaginary part, we obtain at k = 



lmx~\k = 0,lu > 0) = -e 



2 %(n + 2) 1 



d 3 qi d 3 q 2 



(n + 8) 2 SlJ (27r) 3 (27r) 3 
3[n(\qi + g 2 |)(l + n(qi) + n(q 2 )) - n(gi)n(g 2 



4|gi + <?2|gi<?2 

1 + (3/2)(n( 9l ) + n(q 2 ))(l + n(\ gi + q 2 )) 

4|gi + ?2|?i?2 



5(u; + \qi + g 2 | - gi - g 2 ) 



8(u - |gi + g 2 | - gi - g 2 ) 



a; 



> VeT (2.33) 



with Imx _1 (A;, to < 0) = — Irnx -1 (A;, —uj), and where n(x) = l/(e x l T — 1) is the Bose function. 
A similar result can also be obtained for k ^ but we will refrain from displaying it; we will 
limit ourselves to analyzing the simpler k = result. The angular integrals in (|2.33|) can be 
performed and we obtain then 



mix = -e- 



,2n(n + 2) 
(n + 8) 2 



pOO poo 

3 / dqt dq 2 [n(q 1 +q 2 -u)(l + n(q x ) + n(q 2 )) - n(gi)n(g 2 )] 

Jui/2 Ju)2 

1 + ~ W?i) + n(92))(l + n(gi + g 2 - w)) oo > v^T (2.34) 

Somewhat unexpectedly, all of the integrals in ( |2.34| ) can also be performed analytically; 
after a lengthy, but straightforward, computation we obtained a final result which had a 
surprisingly simple form — we got 



ru/2 


nu>2 




+ / 


'qi / c 




Jo 


JW2 — qi 





Imx 1 (0,v) 



,27r(n + 2) 



a; 



(n + 8) 2 

where Li 2 (x) is the dilogarithm function 



n 2 T 2 + 6T 2 U 2 (e~^ 2T ) 
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Li 2 (x) = - f ^ln(l-y). 
Jo y 

For large we have from fl2.35|) 



Imx 1 (0, \u\ — ► oo) 



sgn(u;) 



^2.35) 



(2.36) 



(2.37) 



where 77 = e 2 (n + 2)/ (2(n + 8) 2 ) is the field anomalous dimension; this is precisely the result 
expected || at this order from scaling. For small 00, the result (|2.35 ) taken at face value 
gives us 



Mn + 2) 



[n 



27r 2 T 2 sgn(a;) + 3luT (\u(\lu\/2T) - 1) 



[2.38) 



This singular behavior at small uo is clearly an artifact of taking ( |2.35| ) beyond its regime 
of validity; we expect instead that lmx _1 (0, uo) ~ 10 for small u, but have no direct method 
here for estimating its coefficient. 
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One measure of the strength of the dissipation computed above is the value of 
lmx _1 (0,co> = T), where our expansion is expected to be reliable. This is characterized 
by the damping rate Tr? defined in ( II . 14| ) . From ( |2.35| ) and ( |2.19| ) we have to leading order 
in e 



T 



RT 



13.770249- 



+ 7r 2 + 6Li 2 (e 



-1/2N 



n + 8 



(2.39) 



It is interesting to compare this value with the exact result for the one-dimensional 
transverse- field Ising (n — 1, e = 2), for which we get 



RT 



T 



RT 



2tt 2 



' r(i/i6) \ 
r(i5/i6); 



r 



15 
16 



x 

47T, 

exact value for n 



3.06... 



sin(7r/8) sinh(l/2) 

1,6 

e expansion (2.39) at n = 1, e = 2 



= 2.560527... 
2 



(2.40) 



(the T functions on the right hand side, should not be confused with the damping rate Trt). 
The agreement is quite reasonable, even for e as large as 2. 

It is interesting to compare the ratio of relaxation rates at uj = T (Tr, defined in ( |1.12| )) 
with that at uj = (r^) for the n = 1, d — 1 case, where we have results for both. We 
obtain 



RT 



7T 



X15/16)V 



2sin 2 (7r/16)sinh(l/2) ^(1/16) J 
0.981700183338266.. 



15 
16 



i 

An 



(2.41) 



Notice that the two rates are almost exactly equal, as had been conjectured earlier for d = 2, 
n = 3 in Ref 0. We suspect that the near equality is quite general, and so Trt is always a 
good estimate for Ft in the high T limit (region II of Fig p. 



III. MODELS ABOVE THEIR UPPER CRITICAL DIMENSION 

The computations now follow the same basic strategy as that used for systems below their 
upper critical dimension in Section |T|. The main difference is that the expansion is now in 
terms of the bare value of the irrelevant non-linearity uq, rather than its universal fixed-point 
value. Further, there are no non-trivial renormalizations, and the renormalization constants 
Z, Zi, Z2, Z4 can all be set equal to unity. The results now will have some explicit, non- 
universal, cutoff dependence which cannot be removed by a simple renormalization: this is 
because the T = quantum critical point is above its upper-critical dimension, and the field 
theory is therefore non-renormalizable. 

We will analyze a class of models with the general effective action, S g , of the form 

$1 = |E / tJ^I^'^)! 2 [%) + k 2 + to] + ^ fJ T dr J rfV(,,r), (3.1) 
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where we have the usual Fourier-transformed field 



(f)(k,uj n ) = [ 1/T dr I d d x${x,r)e i(kx -^ T \ (3.2) 



o 



Different choices for M(q, u> n ) describe a variety of physical situations: 

(a) M(u) n ) = u 2 : This obviously corresponds to the action S of a quantum rotor {n > 2) 
or transverse-field Ising (n = 1) model which we have already studied in Section |TJ. It has 
dynamic exponent z = 1 and upper-critical dimension d = 3. 

(b) M{u n ) = —iDuj n : Now S g describes a dilute bose gas with D a constant (analogous 
to the velocity c for S) related to the mass of the bosons. The dynamic critical exponent is 
z — 2, and the upper-critical dimension is d = 2. 

(c) M(uj n ) = D\u n \: In this case, S g describes spin fluctuations in the vicinity of the onset of 
spin-density wave order in a Fermi liquid [l]. Unlike the cases I and II, the T = dynamic 
susceptibility now does not have a quasiparticle pole in the paramagnetic phases, but instead 
has a cut describing the particle-hole continuum. The dynamic critical exponent is z = 2 
and the upper-critical dimension is d — 2. Finally, it must be noted that in this case S g is 



only applicable in the paramagnetic (or Fermi liquid) phase ; a separate action is needed 
within the magnetically ordered phase. 

Note that all of the above choices for M share the property M(0) = 0. Also in all three 
cases the correlation length exponent v = 1/2, and the coupling uq is irrelevant at the 
uq = quantum critical point with scaling dimension — 9 U ; for the models above we have 
9 U — d + z — 4, a relationship which is not always valid. 

Another model of interest is the quantum-critical point describing the onset of ferro- 



magnetism in a Fermi liquid [|TJ. It has recently been pointed out |33[ the effective action 
now contains non-analytic dependencies on the momentum, k, (~ A; d_1 in clean systems and 
~ k d ~~ 2 in random systems) which are present only at T = 0. This singular behavior is 
possible because gapless fermion modes are being integrated out. It is now clearly necessary 
to also account for the T dependence arising from the elimination of the critical fermion 
models. This should be possible using the general methods of this paper, but this issue shall 
not be addressed in this paper. 

Returning to models (a)-(c) above, we perform a perturbation theory in u as described 
in Sections and [TT]. The generalization of ( |2.2| ) to linear order in w is now 



;n + 2\ r d d q ^ 1 T f de 1 . 

^-'o + M— J./ j2^{ T ^ M(e n ) + qi + t + 7 2 ~J^M(e) + q^- ^ 



The susceptibility, defined in (|1.4j) is obtained from R by generalizing ( [2.19| ) to 

= * 2 + %) + R ~ u T (It±2) pigp fl (^ (3.4) 

We have assumed above, and will continue to assume below that 2 < d < 4. The correlation 
length, £, is given at this order by £ 2 = x(0, 0), as in Section |II A 2| . To apply the analog of 
( j2.ip , it is convenient to separate R into the following form: 

in 4- 2\ 

R = t + u f-^-J [Ri + R2 + Rs) (3.5) 
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with 



Ri = 

R 2 = -T 




de 



(2n) d \ ^M(e n ) + q* + t J 2vr M(e) + g 2 + *„, 
d d q ( 1 1 \ 



(2vr) d \q 2 + t q 



2 



^ / (2n) d I 2tt ( M(e) + g 2 + t M (e) + g 2 ) ^'^ 

The integral in R\ is ultraviolet convergent in all d in models (a) and (b), and is convergent 
for d < 4 in model (c). The integral in R 2 is ultraviolet convergent for d < 4. All of the 
ultraviolet divergences have been isolated in R%. For all models (a)-(c) this divergence can 
be separated by a single subtraction which is a linear function to] we write R3 as 




1 1 tn 



M(e) + g 2 + 1 M (e) + g 2 (M(e) + g 2 ) 2 , 
A d d q f de t 



(27r) d J 2vr (M(e) + g 2 )^ 



(3.7) 



The last integral is a cutoff (A) dependent term which has the simplifying feature of being a 
linear (and therefore analytic) function of to- The first integral is now ultraviolet convergent, 
but is a more complicated function of to- In other models additional subtractions involving 
higher powers of to may be necessary at this stage. 

Carrying out all the frequency integrals and summations in Q3.6p and ( |3.7| ) for (along 
with the momentum integration of the last term in ( |3.7| ), we find that in all three models R 
takes the form 

R = t (l - Mo ciA e ") + u c 2 T( 1+e ^T (c 3 ^;) (3.8) 

where c\, c 2 , and C3 are constants, and T is a universal scaling function given by Ri + R 2 + R^ 
but with the last term in ( p.7|) omitted. By examining the limiting behavior of (|3.4j ) and 
( p.8[ ) we can delineate the different physical regimes as shown in Fig [3] P, |l0| , |i"9 . 
(I) T < (t /u ) zu ^ 1+e ^: low T limit of CQFT; paramagnetic phase 

In this regime we are dominated by the t > ground state, with given £ given to leading 
order by its T = value. The subleading temperature dependent corrections are however 
different depending upon whether the argument of the scaling function T(y) is small or 
large. These sub-regimes are therefore: 
(la) T < tg" 

The nature of the corrections depends upon the behavior of T(y — ► 00), which can vary 
considerably from model to model. In models (a) and (b) the t > ground state has a gap, 
so the leading correction will be exponentially small in temperature. Model (c) is a Fermi 
liquid for to > 0, and has power-law corrections in T which will be described in more detail 
below 

(lb) % v < T < f u/{1+euU) 

Now the T-dependent corrections involve T(0), which is always a pure number; so we have 
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£- 2 (T) = C 2 (T = 0) + M c 2 T (1+e " ly)/ ^T(0) + . . . 



(3.9) 



(II) T > |* |W(i+9*«): high T limit of CQFT 

Now T is the most important energy scale and all to-dependent corrections can be neglected. 
The correlation length still obeys (3.S) but with the second term now being larger. 

The transition to the ordered state occurs, as before, at R = 0. To leading order in u , 
T c is given by (p|) to be T c = (\t \ / \u c 2 T (0)) W /W . 

In the remaining presentation we specialize to model (c); the properties of models (a)-(b) 
will be quite similar. The upper critical dimension is d — 2, and we assume we are above it. 
For this case, the explicit result for R is 



R = tn 1 



u A 



d-2 



D 



n 



6 



S, 



ir(d 



+ ^{TD) d ' 2 



n + 2\ ( t 



T 



\DT 



(3.10) 



where T is initially obtained as 

d d q 



T(y) = - [ 

7T J 



(27T) 




2tt 



q +y 



'q 2 + y' 

2tt . 



+ 



n + y 

q 2 



(3.11) 



where if) is the digamma function; it can be verified that the integral over q in (|3.11|) is 
convergent for 2 < d < 4. Now we use the identity if>(z + 1) = if){z) + 1/z to simplify Q3.11|) 
to 



T(l/) 



1 



7T 



(27T) d 




(3.12) 



In this form, it is manifestly clear that T(y) is analytic as a function of y at y = 0, and so 
from ([3.10|) , i? is analytic at t — 0. Indeed the first singularity of (|3.12|) is at y = —2n, and 
( P-12j ) can be used for all y > —2n. This allows us to access the region with t < 0, but 
T > T c (t ). The singularity at y — —2n is of no physical consequence, as it is well within 
the ordered phase. Recall that a similar phenomenon occurred in our earlier analysis of S 
for d < 3 in Section [11 A 1| . 

By evaluating the large y behavior of T(y), we can determine from (3.4) the T-dependent 
corrections in regime la. We find 



T 2 (T) = C 2 (T = 0) + 



u D (w + 2)r((4-rf)/2) 

t (4-d)/2 36(47^/2 



T 



(3.13) 



As noted earlier, the T-dependent correction is a power-law characteristic of a Fermi liquid. 



IV. CONCLUSIONS 



This paper has provided a general strategy computation of finite temperature universal 
crossover functions near quantum-critical points. The strategy can be broken down into 
steps, each step containing distinct physical effects; this separation is an important advan- 
tage of our method. The steps are: 
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(i) Renormalize the T = CQFT to obtain a well-defined quantum theory whose ground 
state, excited states, and scattering amplitudes between them, are known. In principle, this 
information completely specifies the non-zero T properties, and no further renormalizations 
should be necessary. 

(m) Use the information in (i) to integrate out all degrees of freedom with a finite Mat- 
subara frequency, to derive effective action (which could be quite complicated) for the zero 
frequency mode. 

(in) Analyze the effective action by an appropriate technique of classical statistical mechan- 
ics. 

We have applied this method in this paper to a relativistic n-component 4 theory below its 
upper-critical dimension, and to a class of models above their upper-critical dimension. 

We will now comment on the relationship of our results to some earlier work. 

The early work on quantum-critical points J[|,|2j|| studied only the quantum-to-classical 
crossover in the shaded region of Figs |I] and The crossovers in the remainder of the phase 
diagram, and their universal properties, were missed. 

Rasolt et. al. studied the quantum-to-classical crossovers in the dilute Bose gas in d — 3. 
In the present language, these are the crossovers near the T = quantum-critical point at 
chemical potential fi = 0. They described the physics in terms of the Gaussian-Heisenberg 
crossover of 4 field theory in ci-dimensions. This is closely related, but not identical, to 
our description in terms of tricritical functions, as the latter are the universal limit of the 
former when uq <C A A ~ d [32], where A is a momentum cut-off. Our approach properly 
identifies all the non-universal cut-off dependence as due to the T = quantum theory (in 
i?3 in Eqn ( |3.6| )), and shows that the finite temperature corrections are universal (the scaling 
function T in (|3.8|) ). Below the upper-critical dimension, there are no non-universal cut-off 
dependencies, and the use of tricritical crossovers is essential. 

Another popular approach to the study of finite temperature crossovers has been the 
momentum-shell renormalization group (RG), in which the RG equations are T-dependent 
and T is itself scale dependent |1I|,[7|,[8|, [l0|J12]| . This method has been quite useful in identifying 
qualitative features of the crossovers in static and thermodynamic quantities. However, 
quantitative crossover functions have been quite difficult to obtain. We believe this is not 
merely a technical difficulty, but an intrinsic problem with the physical basis of this approach. 
The dynamic consequences of quantum and thermal fluctuations are physically quite distinct, 
and it appears quite unsound to interpolate between them by defining a scale-dependent 
temperature. It is quite clear that such a method will not correctly describe the thermal 
dissipative dynamics. Instead, our point of view is that the RG flows are more properly 
considered as properties of the T = theory, and allow one to define its eigenstates and S 
matrices. The finite T physics is then completely determined by these properties. 

O'Connor and Stephens |22| have used an idea similar to that in the momentum- shell RG, 
but in the framework of the field-theoretic RG. They achieve this by defining some unusual 
renormalization conditions which seem designed to yield (3 functions which are temperature 
dependent. The physical meaning or mathematical justification of these renormalization 
conditions is not clear to us, and our critiques in the previous paragraph apply here too. 
We also note that their results are not systematic expansion in some control parameter 
(many of the terms contain e to all orders), and are not naturally expressed in the terms of 
renormalized T = energy scales which expose the full universality of the physics. 
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Large n expansions have been also been used to study finite temperature properties 
of quantum critical points H9,|34"|. They have the advantage of being uniformly valid over 



the entire phase diagram. Their most extensive application has been in d = 2 ||, where 
T c (t) = for large n. However T c is non-zero for d > 2, and a large n computation then 
gives results consistent with those of this paper. 

We now turn to discussing some open problems and directions for future research. 

The major gap in existing results is a quantitative and systematic theory for the low- 
frequency dynamics. This is an experimentally important question, as the damping rate 
directly determines NMR relaxation rates in two-dimensional quantum antiferromagnets |5j 
Our present approach fails for u < y/eT, but it is possible that a systematic analysis of a 
self-consistent approach, with e a control parameter, can be performed. 

The present paper has avoided discussion of logarithmic corrections in special dimensions, 
either due to d being the upper-critical dimension of the quantum-critical point, or because 
of the logarithmic corrections that appear in weakly-coupled classical theories in d = 2. We 
made this choice to streamline our discussion, but it should not be too difficult to extend 
our results to include these cases. 

Finally, we have already noted that there should be interesting finite T crossovers in 
nearly ferromagnetic Fermi liquids, as some novel non-analytic q dependencies in the effective 



action have recently been pointed out [33 
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APPENDIX A: CLASSICAL TRICRITICAL CROSSOVER FUNCTIONS 

In this section we will tabulate results on the classical tricritical crossover functions 
needed in the region |T — T c (t)\ <C T c (t). We will confine our attention to the crossover 
function, \V c (q,v), appearing in (|1.24j ), for the static susceptibility. In the weak coupling 
region, v <C 1, we can easily expand in a power series in integer powers of v: 

T i/ x 2 fn + 2\ 2r((4-d)/2) o. /A , 

In this region the tricritical crossovers connect smoothly with our e expansion results in the 
region outside |T — T c (t)\ <C T c (t); the equivalent of ( [All ) was already used in (|2.19| ) and 



Within \T — T c (t)\ -C T c (t), v becomes large, and alternative perturbative expansions are 
needed to obtain tricritical crossovers. We will discuss two such methods here. 

The first method is the expansion in 4 — d. The reader may be bothered by our simul- 
taneous use of an expansion in e = 3 — d in the analysis of S in the main part of the paper. 
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However, the two expansions occur in separate calculations and compute entirely different 
crossover functions. They are combined only in the final result, in which the results of one 
appear as arguments of the other. So there is no inconsistency, and the procedure is entirely 
systematic. The result for can be read off from earlier results ||32| ; to leading order in 
4 — d we have 



q 2 + 



1 + 



(n + 8)v 
48tt 2 (4 - d) 



-(n+2)/(ra+8) 



(A2) 



When combined with ( |1.24|) , we see that the critical point is at R — 0; this will change at 
higher orders in 4 — d, when we expect a critical value R ~ [TU / K d ^ 2 ] 2 ^ 4 ~ d \ 

The function ^>c can also be obtained in a large n expansion, with d now arbitrary. 
Taking the large n limit of (|1.23j ) while keeping nv fixed, a straightforward calculation gives 
to leading order 



^c\q,v) = q 2 + Tl(v), 
where 11(f) is determined by the solution of the non-linear equation 

r((4-d)/2) rTT ^ l( d-2)/ 2 



(A3) 



II(f ) + nv 



3(rf-2)(47r) d / 2 



poor 



Finally, it is easily checked that ( |A1| ), ( |A2{ ) and 
mutually overlapping regimes of validity. 



1. (A4) 
all agree with each other in their 



APPENDIX B: COMPUTATIONS FOR S AT T = 



We consider properties of the model S (Eqn ( |1.1| )) at T = and d < 3, in an expansion in 
powers of e = 3 — d. We will compute the renormalized T = parameters which characterize 
the ground state, and appear as arguments of the quantum-critical scaling functions. The 



computations are standard |2^j2J|, and we will be quite brief. 

The renormalization constants Z, Z 4 (and Z 2 ) to the needed order in g in the minimal 
subtraction scheme are: 

n + 2 9 

Z = 1 —g 2 

144e 

n + 2 
Z 2 = l + -—g 
6e 

, n + 8 
Z 4 = l + ——g 
6e 

The fixed point on the /3-function is at g = g* with 



(Bl) 



6e 



n + 



1 + e 



20 + 2n - n 2 
2(n + 8) 2 



(B2) 



We consider the cases t > and t < separately. 
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1. t > 



At T = 0, all properties are "relativistically" invariant, and are most conveniently ex- 
pressed in terms of a Euclidean momentum p = (cj,/c). The renormalized susceptibility x 
takes the form 

X - l (p)=p 2 + t-Z(p 2 ) (B3) 

where £ is the self energy. The quasi-particle pole occurs at p 2 = — A^ which is the solution 
of A 2 — t — £(— A+). The residue at this pole, .4. is given by 



d£ 



(B4) 



P 2 = -A2 

To leading order in g, we can now easily obtain by the usual methods 

a9 / n + 2 \ n + 2 u e r d d+1 p ( 1 1 \ . 

Evaluating this at g = g* we obtain 

A 2 , = /. W) 2 ' (B6) 

where z/ = 1/2 + e(n + 2)/(4(n + 8)) is the correlation length exponent, and there is no 
correction to the prefactor at order e. 

To obtain the leading contribution to A we have to go to order g 2 , where we obtain 



n + 2 2 n + 2 ( fx e g\ 2 d r d d+1 Pl d d+1 
+ 1aa7 9 + ~L$r{~S~^~J dp^J {27TY+ 1 (2tt)< 



_P_2_ 

c2+l 



(Pi + A 2 +) (p| + A^)((p + Pl + p 2 ) 2 + A 2 + ) 



(B7) 

9 2 = -A 2 , 



The integral can be performed by transforming to the usual parametric representation, which 
yields 

n + 2 2 n + 2 2 ^ u \ 26 T 2 (2 - e/2)r(e) ^ ^ , (1 - y)y-^/ i (l - xf/ 2 
A=1+ U4^ 9 -^T 9 [a-J 4 Jo dX L dV (l-x(l-x)(l-y)Y 

(B8) 

Evaluating the integral as a power series in e, we find that the poles in e cancel. Finally, 
replacing g — > g*, we find 

^=(^)"(l + 0.2823615146^*) (B9, 

where the exponent n = (n + 2)e 2 / (2(n + 8) 2 ). 



26 



2. t < 



First we determine the value of Nq = (0). Ordinary bare perturbation theory gives 

d d+1 p ( 1 



'6j*o| 
u 



Up 
4|tnl 



(2tt 



\d+l 



,P 2 + 2|t | 



1 

p 2 



(BIO) 



Re-expressing in terms of the renormalized t and g, and evaluating at g = g*, we find (we 
can ignore the wavefunction renormalization Z at this order) 



o = A* 



+ 8 2\t\ 



2e 



(Bll) 



where the exponent ft = 1/2 — 3e/(2(n + 8)). 

To get the energy scale measuring deviation from criticality, we consider the cases n = 1 
and n > 2 separately: 



a. n = 1 



Bare perturbation theory tells us that xip) ls given by is 

x~\p) =P 2 + t + 



u N 2 u 2 N 2 f d d+1 Pl 



(27T)^((P + Pi) 2 + 2|to|)(^ + 2|t | 



«0 

2 



d d+1 Pi 
(2tt) 



1 



d+l 



v Pi+2|to 



Pi, 



(B12) 



Using ( |B10[ ), re-expressing in terms of the renormalized t and g 1 and evaluating at g = g*, 
we find (again ignoring Z at this order) the energy gap A_ by solving x~ 1 {p 2 = — A 2 ) = 0: 



A^ 



V 2 1 + 



7T 



\/3-3 



2N 
,M 2 



2;/ 



(B13) 



6. n > 2 



In this case we will use the stiffness, p s as a measure of deviation from criticality. We 
compute the transverse susceptibility (measured in a direction orthogonal to the condensate) 
in bare perturbation theory: 



xl\p) =p 2 + t + 



u N 2 



n 2 N 2 
"o iV o 



d d+1 Pl 



1 



(27r) d +i (p + pi) 2 (pl + 2\t \) 
Uq f d d+1 pi ( 1 



+ 



6 J (27r) d +! \p\ + 2\t,\ pi 



(B14) 



Again, we use (P10|) , re-express in terms of the renormalized t and g, and evaluate at g = g*, 
to obtain: 
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x±(p) =p 2 + 



2e\t\ 








n + 8 





log 1 + 



2\t\ 



- 1 



In the region p 2 <C \t\, the above result takes the simple form 



xl\p)=p 2 (l + 



2 n + 8 



We can therefore identify 



2\t\ 



't\V 



2(n + 8) I \/j, 2 



(B15) 



(B16) 



(B17) 



3. Universal ratios 

For completeness, we list here the universal ratios that can be constructed out of the 
t > and t < results. For n — 1, we can take the ratios of the gaps A_ and A + 



A, I 12 



(B18) 



The analog of this ratio for n > 2 is 

Ps 



A'fT 1 V 2e / \" ' 2(n + 8 

A second set of ratios emerges from the ratios of the field scale. Now we have for n = 1 



(B19) 



A^ 2 _ n+_ 
At 1 A ~ 2e 



6 



(B20) 



and for n > 2 



2(n + 8)' 



(B21) 



APPENDIX C: COMPUTATIONS FOR 5 FOR T > 

This appendix will present formal results for the system S (Eqn ( |1 . 1| ) to order at 



nonzero T. In Section [LA] , we outlined how to compute these using a two-step process: (i) 
obtain an effective action <S e // (Eqn ( |1.17|) ) for the u) n = mode; (ii) compute of correlations 
of observables under <S e //. 

First, for future use, let us obtain the value of the mass subtraction m^. Consider the 
susceptibility of the theory in <S e //, but with bare mass ml; to order Uq, this is given in bare 
perturbation theory at T = by 
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d d+1 



Pi 
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P2 
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d d+i„ 
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{2n) d + l {pi + ml) 2 i 
1 



(2n) d + 1 (2n) d + 1 (p\ + ml) {p\ + m§)((p - p x - p 2 ) 2 + m 2 c 



(CI) 



The critical point is determined by the value mo = m§ c at which % 1 (p = 0) = at T = 0. 
Solving ( |CT|) for this condition order by order in uq we obtain 



m, 



0c 



-U 



n + 2\ r d d+1 p 1 9 
+ u 2 
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(27r) d + 1 p 2 



n + 2\ /" d d+1 Pl d d+1 
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(C2) 



In subsequent computations, in all propagators carrying a non-zero frequency, we will insert 
the mass 



m 0c + 1 , 



(C3) 



with ?t?,q c given by (|(J2D , and expand in powers of Uq: all non-zero frequency propagators in 
the resulting expression will therefore have mass t . 

Let us now obtain the couplings in the effective action S e ff for the u n = mode to order 
Mq. By ordinary perturbation theory in the finite frequency modes we obtain 



C 2 (k) = Zk 2 + Zt + u 
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We have implicitly assumed above that Z = 1 + 0(g 2 ) and only written Z where it is need 
for the order g 2 result; we will continue to do this in the remainder of the appendix. In a 
similar manner, we can obtain the value of C^. 



C^k^k^k^-kx -k 2 -k 3 
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where the symbol Sym fc denotes that the expression following it has to be symmetrized 
among the momenta k\, k 2 , fc 3 , fc 4 = — k\ — k 2 — k 3 . All other couplings in S e ff are zero at 
order Uq. 
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We now perform the renormalizations of the super-renormalizable classical theory to 
obtain C 2 and C 4 . First, to order Uq, it is easy to see that C 4 = C4. For C 2 , in addition 
to the tadpole renormalization in ( |1.18| ) there is a two-loop renormalization that has to be 
included for d close to 3 



C 2 (k) = C 2 (k) + 



n + 2\ f d d ki C 4 (k, -k, k x , —hi) 
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(2n) d k\ 
d% d d k 2 C 4 (k,kiJ- 
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-k-kt- k 2 )C 4 (-k, -ki, -k 2 , k + k x + k + 2) 



(2n) d (2n) d (k 2 + T 2 )(A; 2 + T 2 )((k + k x + fc 2 ) 2 + T 2 ) 

This, in fact, completes the set of renormalizations, and there are no new terms that have 
to be accounted for at higher orders in Uq in the super-renormalizable classical theory. To 
avoid an infrared divergence in d — 3, we have performed the two-loop renormalization 
above at an arbitrarily chosen Pauli-Villars mass equal to T. Our results for the coupling 
constants will therefore depend upon this choice of renormalization scheme, but all physical 
observables will be independent of it. Combining ( |U4|) and ( |C6| ) we get 
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The values of the couplings i£, i^, and U now follow directly from (|1 .20 ) and the results 
O, (|C3) above. 

Let us now apply the above procedure to obtain the perturbative result for the dynamic 
susceptibility x(k, iuj n ) at finite T. As discussed in Section [TA] , this result will be valid 
everywhere in the phase diagram of Fig|l|, except in the shaded region \T — T c {t) \ <^ T c (t). 
The simplest way to proceed is to introduce an external source term coupling to the field 
4> and then to proceed in the two-step procedure noted at the beginning of this Appendix. 
We omit the details and state the final result 
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a(e n ) = + RS en ,o + t (l - 5 en>0 ). (CIO) 

Notice that in ( |C9[) we do not expand out the uo dependent expression for R given in ( p. 2D , 
but instead treat R as variable formally independent of uo; this is required by the method 
of Section |l~A| . 

The equations ( |C7| ) and (|C9|), are the main results of this appendix, and will be used in 
the body of the paper. 

In the following subsections of this appendix, we will evaluate the formal results above 
to obtain explicit two- loop expressions for some quantities at t = 0. Our main purpose in 
doing this is to demonstrate the consistency of our approach, by explicitly displaying the 
cancellation of all ultraviolet and infrared divergences and the collapse of the results into 
the scaling forms of Sec. |. 



1. Evaluation of R 



As noted above, we will restrict our results to the critical coupling t = 0. We begin with 
( P?D and the definition ( |1.20| ). Explicitly evaluating out the one-loop contributions in terms 
of the functions introduced in Section [II A 1| (and identity ( |2.16|) ), and expressing in terms 
of the coupling g using (|2.3j), we find to order g 2 
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Notice that the above expression has poles in e multiplying the thermal function -F^(O). 
Consistency requires that these poles must cancel divergences coming out of the two-loop 
frequency summation left unevaluated in ( |U11 ). This is indeed what happens. We can see 
this by adding and subtracting the following expression to ( |C11|) : 
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We absorb the left-hand-side of ( |U12 ) into the unevaluated integrals in ( |CJ1 1| ) . The right- 
hand-side of (|C12| ) has poles in e which precisely cancel the poles in ( |C11|) . Setting g = g* 
(Eqn (|B2"P) and expanding in powers of e, one finds that the fi dependence of ( |C11| ) also 
disappears; in this manner we obtain 
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where the number / arises from the frequency summations in 
evaluating these summations we find 
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where n(q) = l/(e q — 1) is the Bose function at unit temperature. It can be checked by a 
straightforward asymptotic analysis that the combined integrals in ( |C14|) are free of both 
ultraviolet and infrared divergences; we evaluated the integrals numerically and found 
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We also quote the values of the other constants in ( |(J13| ): 
(7(0) = 2.45380858207. . . 
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where 7 = 0.577216 ... is Euler's constant, and £(s) is the Reimann zeta function. 
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2. Evaluation of x(0, 0) 
We can obtain an expression for the static susceptibility, x(0,0), at t — directly from 
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Expressing in terms of g using (2.3), and rearranging terms a bit, we obtain 
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Let us now evaluate some of the integrals in flC18|) to the needed accuracy in e. First, 
we have 



1 



d d qi 



R 



S d+1 J (2n) d q 2 (q 2 + R) 
d d qi 1 



n a/9 ( 1 - 21n2 
2nR {1 - e)/2 1 + e 



S d+1 J (2nY (q 2 + Rf 
1 f d% v 1 

s d+1 J (27r) d ^ (n 2 n + q 2 ) 2 



T 



-1-e 



--G'(O) . 



(C21) 



where the last equation is related to ( 2.16|) . The integrals over momenta in ( p!9|) can be 
performed exactly in d = 3 (which is all we need) and give 
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Now notice that R ~ eT 2 <C T 2 . In this limit we can get the leading result for I 2 simply by 
expanding ( |C22| ) to leading order in R; this leads to 
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The integral in ( |C20| ) can also be evaluated and we find 

J 3 = 2vr 2 ln( J R/T 2 ). 
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We are now ready to assemble all these results into ( |C18| ). Expanding ( |C18| ) in powers of 
g to order g 2 one finds, as expected, that all the poles in e cancel. Setting g = g* (Eqn (|B2| )) 
and expanding in powers of e (while keeping R fixed), one finds that all the \i dependence 
disappears and the resulting expression takes the form 
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We have retained all terms, which after inserting ( |C13|) , are of order e 5 ^ 2 or larger. The three 
loop corrections, of order Mq, contain a contribution like e 3 T/\/R ~ e 5 / 2 , so ( |C25| ) does not 
contain all terms of order e 5//2 ; it does, however, include all terms of order e 2 lne or smaller. 
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3. Evaluation of dx^(k, 0)/dk 2 \ k=0 
we have at t — 0, 
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where a(e n ) is defined in (|C10| ) with t = 0. Now add and subtract the following integral 
from the above 
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with the constant I5 = 0.22657560322 . . .. The poles in e then cancel, and the remainder of 
( C26 ) can be evaluated at e = and g = g*. Let us now define the momentum integral 
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where in the second equation we have transformed to the usual parametric representation. 
Then we can write ( |C26|) in the form 
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It is now not difficult to show that the combination of the summation and integration within 
the square brackets in ( |C29| ) is free of both ultraviolet and infrared divergences. In fact, this 
combination is a dimensionless quantity which is a function only of the dimensionless ratio 
R/T 2 . Now we know from ( |<J13| ) that R/T 2 <ti 1, and in this limit, the term in the square 
bracket in ( |C29| ) is dominated by the single term in the summation with e n = Q n = 0; we 
have therefore 
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FIG. 1. Finite temperature (T) phase diagram, as a function of the tuning parameter t, 
of crossovers near a quantum critical point (t = 0, T = 0) below its upper critical dimension, 
for the case where long-range-order survives at non-zero T. The dashed lines indicate smooth 
crossovers, while the full line is the locus of finite temperature transitions T = T c (t). All crossover 
and phase-transition boundaries scale as T ~ \t\ zu , where z is the dynamic exponent and v is 
the correlation length exponent. For the action S (Eqn (|l.l| )) this paper provides a systematic 
expansion of observables with e = 3—d the control parameter. In region II the expansion is in powers 
of yfe (e = 3—d), with additional factors of lne; in the shaded region (defined by \T—T c (t)\ <C T c (t)) 
we provide an e-expansion of coupling constants, which then become arguments of previously known 
classical, tricritical, crossover functions. Elsewhere, all observables can be obtained in an expansion 
in integer powers of e. For T / 0, all observables are analytic as a function of t at t = 0; in the 
region t < 0, T > T c (t), our results are obtained by analytic continuation from the t > results. All 
properties of the phase diagram are described by the continuum quantum field theory associated 
with the quantum critical point, which is in its low T limit in regions I and III, and in its high T 
limit in region II. The condition determining boundary of region III is similar to that determining 
T c to within constants of order unity that we are free to choose, and we have chosen region III to 



extend to both sides of T = T c (t). For more discussion on the regions see Section II A 2 
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FIG. 3. As in Fig |l], but for the case where the quantum critical point is above its upper-critical 
dimension. Now properties are controlled by the value of the coupling uo, associated with the least 
irrelevant operator involving quantum-mechanical interactions; we denote its scaling dimension as 
On (0 U > 0). The crossover boundaries on either side of region II, and the phase transition line 
T = T c (t), now scale as T ~ (\t\/uo) zu ^ 1+9uU \ The boundary between regions la and lb scales as 
T ~ t zv . The expansion is now in powers of the bare coupling uq, and as before, classical tricritical 
crossover functions are needed in the shaded region \T — T c (t)\ <C T c (t). For more discussion on 
the regions see the discussion below Eqn ( |3.8| ). 
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FIG. 4. A plot of the universal scaling function T(y), denned in ( |3.12| ), for d = 3. Notice that 
it is smooth at y = 0. It is analytic for all real y > —2ir. 
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